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I. INTRODUCTION 



A. BACKGROUND 

In 1998, Sandia National Laboratories (Sandia) contracted the 
Naval Postgraduate School (NPS) to study the pressure field associated 
with a proposed aerodynamic fairing, designed to fit beneath the forward 
fuselage section of the Altus II Unmanned Air Vehicle (UAV). The 
purpose for this analysis arose from the requirement to fit an Ultraviolet 
(UV) laser remote sensing payload into the forward fuselage section of 
the Altus II UAV. The dimensions of the UV laser telescope necessitated 
that a portion of the telescope protrude into the freestream air beneath the 
originally configured Altus II airframe. To reduce drag and prevent 
unacceptable levels of telescope vibration, two proposals were developed 
to design an aerodynamic fairing that would facilitate the installation of 
this UV telescope. See Figure 1 for an illustration of the Altus II. 




Figure 1 . Altus II UAV 
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The first proposal was to design a fairing with an open cavity 
window, through which the UV telescope would view. This idea had the 
advantage of minimizing the optical distortion associated with having the 
telescope view through an optical window. One disadvantage of the open 
cavity fairing design was that vibrations associated with the open cavity 
flowfield could generate considerable stress on the telescope. [Ref. 1] 

This stress would in turn impose nearly continuous strain on the 
telescope’s gimbal system designed to aim and hold it onto the desired 
target. 

The second proposal was to design a fairing with a transparent 
optical window, through which the UV telescope would view. The 
engineers at Sandia concluded that a high quality optical window could be 
designed to sufficiently minimize the distortions in the UV telescope 
imagery. One advantage of the closed fairing design was that the strains 
imposed on the telescope’s aiming/holding system would be greatly 




Figure 2. Proposed Aerodynamic Fairing 
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reduced. When all of these considerations were taken into account, the 
decision was made to proceed forward with the closed fairing design. 
Figure 2, provided by Valley Engineering Group (VEG), illustrates the 
proposed shape of the fairing. 

B. PURPOSE 

The original purpose of the fairing analysis was to measure the 
open cavity flow field forces exerted on the proposed gimbaled system 
used to align and hold the telescope in position. The forces and 
vibrations imposed on the telescope in the open cavity design would have 
severely restricted the design of the telescope’s aiming/holding system. 
After the decision was made by Sandia to change to the second proposal 
for the fairing design, the purpose of analysis changed. 

Using the closed fairing with an optical window as the baseline 
design, the purpose of this analysis was to estimate the pressure-induced 
forces on the optical window and to estimate the pressure characteristics 
along the longitudinal centerline axis of the Altus II ’s upper fuselage. 
Knowing the approximate force exerted on the window will enable design 
engineers to fabricate a window possessing minimum weight with 
acceptable distortion and optical characteristics. 
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II. TEST EQUIPMENT DESCRIPTION 



A. THE MODEL 

1. Model Construction 

The wind tunnel model used for this analysis was a one-half scale 
likeness of the forward fuselage section of the Altus II UAV, with the 
proposed fairing attached. The foundation of the model was cut from a 
high-quality pine plywood board upon which the upper fuselage and lower 
fairing surfaces were mounted. To form the shape of these upper fuselage 
and proposed fairing surfaces as closely as possible to the graphic 
depiction provided by VEG, cross-sectional ribs were affixed to the 
baseboard to aid in developing the final surface shapes (Figure 3). 




Figure 3. Fairing Shape Before Fiberglass Covering 
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Blue polystyrene foam, commonly used in fiberglass hand lay-up 
construction, was inserted between each of the fairing’s cross-sectional 
ribs. The blue foam was then sanded down to the level of the cross- 
sectional ribs with the contours in the inter-rib regions ‘interpolated’ by 
visual observation. Once the blue foam shapes were finished, three layers 
of fiberglass were applied to the foam core. After the fiberglass set, the 
foam inside the fairing was removed to facilitate the installation of the 
pressure ports and the associated Tygon tubing. See Figure 4 for the 
finished fairing. 




Figure 4. Painted Fiberglass Fairing 



2. Model Instrumentation 

The principal area of interest on the aerodynamic fairing was that 
area where the optical window was to be mounted. Consequently, the 
entire window surface was instrumented with a pattern of pressure ports. 
To reduce the total number of pressure ports on the fairing, only the right 
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side of the forward and aft sections was instrumented for pressure 
measurements. The pressure ports were divided into three sections: the 
forward section (A Section - 40 ports), the middle section (B Section - 45 
ports) and the aft section (C Section - 40 ports). See Figure 5 for details 
regarding the model port numbers. Dividing the pressure port pattern into 
these three sections was necessary because the Scanivalve’s pressure 
manifold could only accommodate 48 pressure measurements at one time. 



Pressure Port Locations 




Figure 5. Pressure Port Diagram 

Once the holes were drilled for each pressure port, brass fittings 
were individually glued into every port. Each brass fitting was tailored so 
that surface flow disruptions on the model were virtually eliminated. 

Inside the fairing, Tygon tubing was attached to each pressure port fitting. 
Initially, the entire Tygon tubing bundle was divided into two smaller 
bundles to enable the tubing to exit the model and wind tunnel. Once the 
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two bundles were outside the tunnel, they were divided again to reform 
three tubing bundles. These three bundles corresponded to the three 
aforementioned pressure port sections on the model: the forward (A), 
middle (B) and aft (C) sections. 

B. THE WIND TUNNEL 

The NPS 3.5’ x 5.0’ Academic Wind Tunnel (AWT) was used to test 
the proposed aerodynamic fairing. The vertically oriented closed circuit 
AWT shown in Figure 6 measures approximately 62.7 feet in length, 33.0 
feet in height and 15.0 feet in width with a 14.4 ft 2 (2079 in 2 ) test section. 
[Ref. 2] The AWT was originally powered by two 150-hp electric motors 
that drove two counter-rotating three-bladed variable-pitch fans. As a 
result of a previous incident involving the fan ingestion of a hand tool, 
the AWT currently uses only one 150-hp electric motor to drive a single 
fan. Single motor operation gives the AWT an approximate maximum test 
speed of 145 KTAS. [Ref. 2] 




Figure 6. 3.5’ x 5.0’ Academic Wind Tunnel 
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The Altus II UAV operates at low subsonic speeds, typically well 
below 120 knots. For this reason, the AWT is well suited for testing a 
model of this type of flight vehicle. The frontal cross-sectional area of 
the Altus UAV model measures approximately 227 square inches, as 
shown in Figure 7. Based on this cross-sectional area, the blockage 
equals 10.9%. The model is relatively large for the given test section; the 
associated blockage effects are discussed in Chapter III Section B. The 
model scale was chosen to minimize Reynolds number effects. 




Figure 7. Frontal View of Model in Test Section 



C. THE DATA COLLECTION EQUIPMENT 

Data collection was accomplished by integrating a Scanivalve 
pressure measurement device with a PC-based National Instruments 
Labview data acquisition program. See Figure 8 for a schematic diagram 
of the test equipment set-up. 
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Figure 8. Schematic of Test Equipment Set-up 



1. The Scanivalve 

The Scanivalve is a differential pressure measurement device that 
used a calibrated transducer to measure the difference between a test 
pressure and a reference pressure (ambient pressure at the Scanivalve’s 
rotary switch). The Scanivalve used in this experiment had a +/- 2.5 psid 
transducer. It mechanically cycled between and measures pressures from 
48 individual ports located on its pressure manifold. Combined with a 
compatible data acquisition system, the Scanivalve was an effective low- 
cost means of obtaining low speed wind tunnel pressure measurements. 

2. The Scanivalve Digital Interface Unit (SDIU) 

The SDIU was an analog-to-digital converter that converted analog 
voltage measurements from the Scanivalve (corresponding to discrete 
differential pressure measurements) into digital voltage readings. It also 
enabled an operator to select individual numbered pressure ports for 
measurement. Additionally, the SDIU is equipped with a GPIB link which 
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enabled it to be operated remotely by an appropriately configured PC. A 
486 PC using NI Labview software controlled operated the SDIU. 

3. The Scanivalve Solenoid Controller 

The Controller provided 28 volt DC electrical power to the 
Scanivalve unit. It also enabled an operator to manually step through or 
reset-to-home the Scanivalve’s pressure ports. This was one of three units 
interfaced together with the Junction Box discussed below. 

4. The Signal Conditioner 

The Signal Conditioner enabled an operator to zero and span the 
SDIU voltage readings. Additionally, it provided the excitation voltage 
for the Scanivalve. To calibrate the SDIU, a water manometer was used to 
calibrate the transducer at 30 cm of H 2 0 differential pressure. 

5. The PC and National Instruments Labview 

The software program used to drive the experimental equipment and 
collect the pressure data was taken from LT Greco’s thesis work. [Ref. 5] 
The NI Labview VI used for his experiment also employed a Scanivalve to 
take pressure measurements. The only modification made to the VI was a 
small addition to the code that enabled the pressure measurement data to 
be stored on a PC 3.5 inch floppy disk. 

6. The Junction Box 

A junction box that interfaced the Scanivalve, the SDIU and the 
Scanivalve Solenoid Controller was fabricated for this experiment. Using 
design drawings provided by Dr. Hobson, this junction box was 
constructed to the specifications of another junction box used in the 
experimental set-up of LT Greco. [Ref. 5] Properly connected to the other 
experimental apparatus, this junction box facilitated the 
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PC/SDIU/Scanivalve interface that enabled the pressure measurement 
system to function as required. 
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III. RESULTS AND DISCUSSION 



A. TEST PROCEDURE AND RESULTS 

After the construction of the wind tunnel model was complete, it 
was mounted inverted in the test section of the AWT. As shown in Figure 
9, the model was supported by two bayonet mounts and one aft-center 
mount, and was free to pivot about its vertical and lateral axes. The two 
bayonet mounts were attached to an internal pivoting mechanism that 
allowed up to five degrees of positive/negative sideslip. The aft-center 
mount was constructed with a screw-post that allowed up to four degrees 
of positive/negative AOA. 




Figure 9. Side View of Model in Test Section 
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lo Tunnel Operation Procedure 

Several preliminary wind tunnel test runs were conducted to ensure 
that the model and support structures were robust enough to handle the 
test velocities. At wind tunnel velocities in excess of AP = 10 cm of H 2 O 
pressure (as measured by the H 2 O manometer at the AWT Control Panel, 
shown in Figure 10), model vibration was observed. AP measured the 
difference in pressure between a Kiel probe and a static ring, located just 
prior to the test section. The conversion from indicated differential 
pressure (AP) at the AWT Control Panel to test section dynamic pressure 
is discussed in Section B of this chapter. 

Two wind tunnel tests were conducted at test section velocities 
corresponding to AP = 5 cm and AP = 10 cm H 2 O to determine whether 
Reynolds number significantly affected the pressure measurements. 




Figure 10. AWT Control Panel and Water Manometer 
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See Table 1 for conversion between water pressure, pound/foot 2 and 
airflow velocities. 



&P(cmU2b) 


Q(psf) 


Q (psi) 


V (feet/sec) 


V (KTAS) 


5 


10.2 


0.071 


93.0 


55.1 


10 


20.5 


0.142 


131.5 


77.9 


15 


30.7 


0.213 


161.0 


95.4 


20 


40.9 


0.284 


186.0 


110.2 


25 


51.2 


0.355 


207.9 


123.2 



Note : density assumed sea-level standard 



Table 1. Some Common Conversions for Low-Speed Wind Tunnel Tests 



See Appendix A for the tunnel operation checklist used to conduct 
each test run. In brief, the test equipment was powered up, the 
SDIU/Scanivalve was calibrated with a stand-alone water manometer, the 
model was set to the desired AOA/Sideslip (a/p) and the tunnel/test 
section was inspected for foreign or loose objects. Next, the wind tunnel 
was powered-up and the flow velocity was steadily increased to AP = 10 
cm of H 2 O pressure. Once the tunnel velocity stabilized, the NI Labview 
program, which measured and recorded the data taken from the first 48 
pressure ports, was enabled. After the initial 48 port measurements were 
complete, the first 48-port manifold was removed and replaced with the 
second 48-port manifold while the tunnel maintained its test velocity. 

This procedure was repeated for each of the three 48-port manifolds until 
the data for all 125 (40 + 45 + 40) model ports was measured and 
recorded. Figure 11 illustrates the three 48-port pressure manifolds 
adjacent to the Scanivalve. 

2. Reynolds Number Effect on C p 

Since it was observed that excessive model vibrations might occur 
at full Reynolds number testing, the effects of reduced Reynolds number 
testing was observed. Two test-runs, one at AP = 5 cm and another at AP 
= 10 cm H 2 O water pressure, were conducted to observe the effects on the 
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Figure 1 1 . The Three 48-Port Pressure Manifolds 




Figure 12. Longitudinal Axis of Fairing for Reynolds Number Effects Tests 
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pressure measurements taken along the longitudinal axis of the model. 
See Figure 12 for the longitudinal measurement axis used for these 
Reynolds number effects tests. 



Cp vs Longitudinal Position 
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Figure 13. Reynolds Number Effects Test - Cp vs Longitudinal Station Position 

Figure 13 illustrates that a reduced Reynolds number has little 
overall effect on the Cp measurements. The pressure coefficient C P is 
given by: 



where p*, and q« are freestream ambient pressure and dynamic pressure, 
respectively. Although some variation in the measurements exists 
between stations four through eight, the difference in these measured C P ’s 
is negligible for this type of analysis. The overall intent of these Cp 
measurements is to estimate the average force exerted on the optical 
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window area (longitudinal stations 7-15 in Figure 13). It was decided to 
run at sub-scale Reynolds number to avoid model aerostructural concerns. 
The Reynolds number for the 18 test cases is determined from 

= where the absolute viscosity of the air is given by Sutherland’s 



equation M = 2.27x1 O' 



(lbf s/ft 2 ). [Ref. 4] The characteristic 



^T + 198.61 

length for the Reynolds number is taken as the model’s fairing length, 



£ = 4 ft. Standard sea-level density p is assumed and the test velocity is 
determined from Q-bar, as discussed in section B sub-section 1 of this 
chapter. Based on these assumptions, test case Reynolds numbers equaled 
3.1 million. Based on a standard 5000-ft altitude flight at 70 KTAS, the 
actual full-scale Reynolds number equals 6.4 million. 



B. WINDOW PRESSURE MEASURMENTS 

The first step in obtaining the total force exerted on the optical 
window was to determine the local pressure coefficient C p at each port 
location on the model. One calibration and two corrections were applied 
to the raw experimental data to obtain the C p ’s. The wind tunnel 
calibration related test section average dynamic pressure Q-bar to 
measured ^P, as displayed on the AWT Control Panel H 2 0 manometer. 
The two corrections took into account test section blockage effects due 
the model size and zero voltage readings from the Scanivalve reference 
pressure. 



1. Data Calibration, Reduction and Corrections 
a. Tunnel Q-bar Calibration 

In Reference 2, the author calibrated the wind tunnel test 
section Q-bar. After conducting a pressure survey of the AWT test 
section, a relation between indicated pressure (a?) and calibrated dynamic 
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pressure (Q-bar) was determined. Figure 14 illustrates the results from 
this calibration. The ratio of the calibrated (Q-bar) to the measured 
(as read from the AWT water manometer) is a linear function. After the 
raw pressure data was collected in the Microsoft Excel spreadsheets, this 
calibration was applied to all of the other pressure measurements. 

Academic Wind Tunnel Calibration 
16 October 1992 




Figure 14. Academic Wind Tunnel Calibration 

All of the test cases were conducted at approximately = 10 
cm H^O (20.5 psf), which after unit conversion and calibration, resulted 
in Q-bar = 20.7 psf. 

b. Test Section Blockage Correction 

The second modification to the raw pressure data involved the 
correction for test section blockage effects. Authors Rae and Pope 
propose a simple correction factor that applies to models of ‘unusual 
shapes’: [Ref. 6] 
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1 Model Frontal Area AV 

£ = £ = 

4 Test Section Area V v 

g is the ratio of the velocity change due to blockage over uncorrected 
velocity. Rae and Pope state that “a maximum ratio of model frontal area 
to test section cross-sectional area of 7.5% should probably be used, 
unless errors of several percent can be accepted.” [Ref 6] The model 
blockage in the current tests is 10.9%. Since average pressure forces are 
the desired result of this analysis, it was decided that correction errors of 
“several percent” are acceptable. 

c. Zero Pressure Differential Correction 

When the pressure measurement data was observed, a small 
inconsistency was noted in the data. Ambient atmospheric pressure was 
not always measured as zero volts. Scanivalve pressure manifold port #1 
was open to the ambient atmosphere, which is also the Scanivalve’s 
reference pressure when taking the other differential pressure 
measurements. In theory, a properly calibrated Scanivalve/SDIU should 
always read zero volts for port #1, since any variation in ambient 
atmospheric pressure would be sensed by both port #1 and the 
Scanivalve’s reference port located near its rotary switch. Since both of 
these ports measure ambient atmospheric pressure, the measured pressure 
difference between them should equal exactly zero. 

During the various wind tunnels runs, the Scanivalve/SDIU 
was calibrated with a stand-alone water manometer for each testing 
session. Each testing session lasted about one hour and measured between 
six and nine Scanvalve pressure manifolds. Three manifold measurements 
were required to measure the entire model in one configuration. During 
calibration, the SDIU was zeroed for ambient pressure and set to 
approximately -1.78 Volts for 30 cm of H 2 0 pressure on the manometer. 
The signal conditioner was used to set this zero and span. Setting the 
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span (of -1.78 Volts) sometimes slightly shifted the zero pressure reading, 
so an iterative process of zero and span setting was used to confirm that 
each reading was stable. 

During the test runs a small measurement drift occurred for 
the ambient atmospheric pressure reading. To correct this problem, the 
voltage (and corresponding pressure) reading of Scanivalve port #1 was 
subtracted from all of the other pressure port measurements for each 
individual test run. This action applied an average correction to the 
measured data, which removed the error induced by the measurement drift 
of the Scanivalve’s reference pressure. 

2. Pressure Coefficient Profiles on the Fairing 

This section contains Figures 15 through 32, which depict the 
contour plots of the measured C p ’s on the aerodynamic fairing. Appendix 
C contains the MATLAB program code used to convert the Microsoft 
Excel spreadsheet data into these pictorial C p representations. Table 2 
presents a summary list of the 18 test cases run at the specified values of 

a and p- 



T est Case # 


"ADA (deg? 


Sideslip (deg) 


lest Case # 


AO A {’deg) 


Sideslip (deg) 


1 


0 


0 


1U 


-4 


-2.5 


2 


-2 


0 


ft 


-4 


-5 


3 


2 


0 


12 


-4 


2.5 


4 


4 


0 


T3 


-4 


5 


5 


-4 


0 


T4 


4 


5 


6 


0 


-2.5 


15 


4 


2.5 


7 


0 


^5 


T6 


4 


-2.5 


8 


0 


Z5 


T7 


4 


-5 


g 


0 


5 


T8 


0 


0 



Note: test case¥T8 conducted with antenna installed 



Table 2. List of a and p Test Cases 



Included in the discussion of each test case is the pressure-induced force 
on the window. These force calculations assume a flight condition of 100 
KTAS at sea level and are discussed in more detail in the next section. 
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Model Dimension, inches Model Dimension, inches 



Pressure Coefficients 




Figure 15. Case 1 * Cp at 0° a and 0° p 



Pressure Coefficients 




Figure 16. Case 2 - Cp at -2° a and 0° p 
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Model Dimension, inches Model Dimension, inches 




Figure 17. Case 3 - Cp at +2° a and 0° (3 



Pressure Coefficients 





Figure 1 8. Case 4 - Cp at +4° a and 0° P 
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Model Dimension, inches Model Dimension, inches 



Pressure Coefficients 




Model Dimension, inches 




0.1 

0 

- 0.1 

- 0.2 

- 0.3 

- 0.4 

- 0.5 



Figure 19. Case 5 - Cp at -4° a and 0° p 



Pressure Coefficients 




Model Dimension, inches 




Figure 20. Case 6 - C P at 0° a and -2.5° p 
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3 3 2- §. 2 Model Dimension, inches 



Pressure Coefficients 





Figure 21. Case 7 - Cp at 0° a and -5° p 



Pressure Coefficients 




Model Dimension, inches 




Figure 22. Case 8 - Cp at 0° a and +2.5° p 
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Model Dimension, inches Model Dimension, inches 



Pressure Coefficients 




Model Dimension, inches 




Figure 23. Case 9 - Cp at 0° a and +5° p 



Pressure Coefficients 





Figure 24. Case 10 - Cp at -4° a and -2.5° P 
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Model Dimension, inches Model Dimension, inches 



Pressure Coefficients 




Model Dimension, inches 



Figure 25. Case 11 - Cp at -4° a and -5° (3 



Pressure Coefficients 




Model Dimension, inches 



Figure 26. Case 12 - Cp at -4° a and +2.5° P 
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Model Dimension, inches Model Dimension, inches 



Pressure Coefficients 

0.1 
0 

- 0.1 
- 0.2 
- 0.3 
- 0.4 
- 0.5 

-30 -20 -10 0 10 20 30 

Model Dimension, inches 




Figure 27. Case 13 - Cp at -4° a and +5° P 



Pressure Coefficients 





- - - 0.2 



Figure 28. Case 14 - Cp at +4° a and +5° p 
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Model Dimension, inches Model Dimension, inches 



Pressure Coefficients 




Model Dimension, inches 



Figure 29. Case 15 - Cp at +4° a and +2.5° p 



Pressure Coefficients 




Figure 30. Case 16 - Cp at +4° a and -2.5° P 
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Model Dimension, inches Model Dimension, inches 



Pressure Coefficients 




Model Dimension, inches 




0.1 

0 

- 0.1 

- 0.2 

- 0.3 

- 0.4 

- 0.5 



Figure 3 1 . Case 1 7 - Cp at +4° a and -5° P 



Pressure Coefficients 




-30 -20 -10 0 10 20 30 

Model Dimension, inches 

Figure 32. Case 1 8 - Cp at 0° a and 0° p with Antenna 



30 



Figure 15 illustrates the results of the baseline test case, which was 
conducted at 0° a and 0° p. The stagnation line (Cp =1) occurred near the 
model leading edge, located well before the pressure measurement area. 
Consequently, the pressure contours do not depict it. The highest surface 
pressures were measured immediately aft of the leading-edge area and the 
lowest surface pressures were located just forward of the window. The 
window area experienced a relatively constant pressure field of Cp~ - 
0.35. Additionally, there were two small areas of local low pressure (Cp ~ 
-0.5) located aft of the window in line with the window longitudinal 
edges. See position (±5, -10) on Figure 14. The total force on the 
window for test case 1, extrapolated to full scale for the nominal flight 
condition described, was 51.3 lbf suction. 

Figure 16 illustrates the next test case, which was conducted at -2° 
a and 0° p. Compared to Figure 15, the leading-edge area of high 
pressure remained largely unchanged. The local suction region forward of 
the window was also similar to the previous case, except that the area aft 
this region experienced a greater increasing pressure gradient across the 
leading edge of the window. After this increasing gradient, the remaining 
window area experienced a relatively constant pressure of Cp ~ -0.3. The 
total force on the window for test case 2 was 40.4 lbf suction. 

Figure 17 illustrates test case 3, which was conducted at +2° a and 
0° p. The leading-edge area of high pressure was approximately twice as 
large as the two previous test cases. The increasing pressure gradient 
associated with the local suction region forward of the window was less 
than the previous two test cases. This case had a Cp ~ -0.4, vice a Cp ~ - 
0.5 from the two previous test cases, for this local suction region. The 
center of the window area experienced a relatively constant pressure of C P 
~ -0.17, compared to -0.35 and —0.30 from test cases 1 and 2, 
respectively. The trailing edge of the window experienced a decreasing 
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pressure gradient from C P ~ -0.2 to C P ~ -0.35. The total force on the 
window for test case 3 was 40.4 lbf suction. 

Figure 18 illustrates test case 4, which was conducted at +4° a and 
0° p. The leading-edge area of high pressure was similar to test case 3. 
The local suction region forward of the window was also similar to test 
case 3, except C P ~ -0.43 instead of C P ~ -0.4. The increasing pressure 
gradient located aft of this suction region was less than test case 3. The 
center window area experienced a relatively constant pressure of C P ~ - 
0.25, compared to -0.17 from test case 3. The total force on the window 
for test case 4 was 37.0 lbf suction. 

Figure 19 illustrates test case 5, which was conducted at -4° a and 
0° p. The leading-edge area of high pressure (C P > 0) was the smallest 
yet, as compared to the previous test cases. The local suction region 
forward of the window was the strongest of the test cases discussed so far 
with a C P < -0.5. The center window area experienced a relatively 
constant pressure of C P ~ -0.35. There was no significant pressure 
gradient at the window trailing edge, as in test case 3. The total force on 
the window for test case 5 was 48.7 lbf suction. 

Figure 20 illustrates test case 6, which was conducted at 0° a and - 
2.5° p. This test case was the first in the series of non-zero sideslip runs. 
The leading-edge area of high pressure is skewed slightly in the direction 
of the sideslip. The model was tested upside down, so negative sideslip 
angles appear as relative wind coming from the right. The local suction 
region located forward of the window is similar to test case 1, but it was 
skewed slightly due to the -2.5° p. The center window area experienced a 
relatively constant pressure of C P ~ -0.27. Located at position (-4, -12) 
on Figure 20, a small local suction region developed due to the -2.5° p. 
The total force on the window for test case 6 was 41.1 lbf suction. 

Figure 21 illustrates test case 7, which was conducted at 0° a and - 
5° p. This test case was similar to the previous one, except the skew 
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angle of the leading-edge area of high pressure was more pronounced. 

The local suction region located forward of the window was also more 
skewed than test case 6. The center window area experienced a relatively 
constant pressure of Cp ~ -0.30. The leeward local suction region at 
position (-4,-12) increased in strength, compared to test case 6. The total 
force on the- window for test case 7 was 44.9 lbf suction. 

Figure 22 illustrates test case 8, which was conducted at 0° a and 
+2.5° (3. This run was the mirror image of test case 6. The results were 
similar to test case 6, except that the skew angle of the pressure field was 
in the opposite direction due to +2.5° p, vice -2.5° p. The total force on 
the window for test case 8 was 46.2 lbf suction. 

Figure 23 illustrates test case 9, which was conducted at 0° a and 
+5° p. This run was the mirror image of test case 7. The results were 
similar to test case 7, except that the skew angle of the pressure field was 
in the opposite direction due to +5° p, vice -5° p. The total force on the 
window for test case 9 was 50.7 lbf suction. 

Figure 24 illustrates test case 10, which was conducted at -4° a and 
-2.5° p. This test case was the first in the series of ones examining the 
combination of non-zero a and p. The leading-edge area of high pressure 
was relatively small (similar to test case 5) and also had the skew angle 
due to the -2.5° p. The local suction region located in front of the 
window was strong (similar to test case 5), but was more skewed than test 
case 6. Apparently, the -4° a in combination with the -2.5° p increased 
the skew angle effect of the pressure field. The center window area 
experienced a relatively constant pressure of Cp ~ -0.34. The total force 
on the window for test case 10 was 52.2 lbf suction. 

Figure 25 illustrates test case 11, which was conducted at -4° a and 
-5° p. At the maximum negative a and p, these results were similar to 
test case 10. The skew angle effect was more pronounced than test case 
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10 due to the -5° (3, vice the -2.5° p. Otherwise, the results of test cases 

10 and 11 were very similar. The total force on the window for test case 

1 1 was 51.1 lbf suction. 

Figure 26 illustrates test case 12, which was conducted at -4° a and 
+2.5° p. This run was the mirror image of test case 10. The results were 
similar to test case 10, except that the skew angle of the pressure field 
was in the opposite direction due to +2.5° P, vice -2.5° p. The total force 
on the window for test case 12 was 50.5 lbf suction. 

Figure 27 illustrates test case 13, which was conducted at -4° a and 
+5° p. This run was the mirror image of test case 11. The results were 
similar to test case 11, except that the skew angle of the pressure field 
was in the opposite direction due to +5° P, vice -5° p. The total force on 
the window for test case 13 was 53.2 lbf suction. 

Figure 28 illustrates test case 14, which was conducted at +4° a and 
+5° p. The leading-edge area of high pressure was similar to test case 4, 
except that it was skewed similar to test case 9. The local suction region 
located in front of the window was similar to test case 4 in strength with 
Cp ~ -0.43. The center window area experienced a relatively constant 
pressure of C P ~ -0.27. The +5’ p caused a local suction region on the 
leeward surface of the fairing located at position (+5,-12) with a pressure 
of C P ~ -0.44. The total force on the window for test case 14 was 45.2 lbf 
suction. 

Figure 29 illustrates test case 15, which was conducted at +4° a and 
+2.5° p. The results of this run were similar to test case 14, except that 
the skew angle was less pronounced because of the +2.5° p, vice +5° p. 
The center window area experienced the same relatively constant pressure 
as test case 14. The total force on the window for test case 15 was 43.8 
lbf suction. 
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Figure 30 illustrates test case 16, which was conducted at +4° a and 
-2.5° p. This run was the mirror image of test case 15. The results were 
similar to test case 15, except that the skew angle of the pressure field 
was in the opposite direction due to -2.5° P, vice +2.5° p. The total force 
on the window for test case 16 was 40.3 lbf suction. 

Figure 31 illustrates test case 17, which was conducted at +4° a and 
-5° p. This run was the mirror image of test case 14. The results were 
similar to test case 14, except that the skew angle of the pressure field 
was in the opposite direction due to -5° p, vice +5° p. The total force on 
the window for test case 17 was 43.2 lbf suction. 

Figure 32 illustrates test case 18, which was conducted at 0° a and - 
0° P with the external antenna installed. The purpose of this run was to 
determine what effect, if any, the presence of the external antenna had on 
the measured surface pressures. When compared to test case 1, which was 
conducted at the same a and P with no external antenna installed, these 
two pressure contour patterns were very similar. There did exist a very 
small asymmetry in the pressure field on the antenna side of the fairing, 
but this had a minimal effect on the window pressures. 

The model antenna was not constructed at the time the first 17 test 
cases were conducted. Given the results of test case 18, the presence of 
the external antenna was deemed as insignificant as it relates to the 
surface pressure measurements on the window. Therefore, the previous 17 
test cases were not repeated with the external antenna installed. The total 
force on the window for test case 18 was 46.9 lbf suction. 

In summary, as a increased, the contour lines where suction 
pressure began (Cp = 0) moved aft toward the optical window. When P 
was non-zero, the asymmetry of the pressure contour patterns on the lee 
corner of the fairing illustrated a local suction region. This agrees with 
what one would intuitively expect. For all 18 test cases, relatively 
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constant pressures were measured across the aft one-third of the fairing. 
This indicates that the flow was probably separated in this aft region of 
the fairing. 

Recall that only the right side of the forward and aft sections of the 
model was instrumented for pressure measurements. Left side data was 
mirrored from the appropriate right side data, as determined by the 
Sideslip angle (3. For P=0, current test case data was mirrored from the 
right to the left side. For non-zero p, data from the opposite-sign P test 
case was mirrored from the right side to the left. 

The optical window area experiences relatively constant Cp for 
every test case configuration. This observation is significant in that it 
suggests that the position of the optical window is relatively well suited, 
in terms of minimizing distortions due to large pressure gradients across 
the surface of the glass. Test cases 1 through 17 were conducted with no 
external antenna installed on the model. Test case 18 was conducted with 
the external antenna installed on the model. It was determined that the 
measured Cp variation due to the presence of the antenna was negligible. 
Therefore, test case configurations 1 through 17 were not repeated with 
the antenna installed. 

3. Total Pressure Induced Forces on the Fairing Window 

To determine a representative total force exerted on the optical 
window during flight, a flight condition of 100 KTAS at sea level was 
assumed. Q-bar at this condition equals 33.9 psf. Given this assumption, 
the calculated Cp at each pressure port on the window was converted to a 
local pressure and integrated across the area of the window. Summing 
these local pressures across each of their respective rectangular elemental 
areas results in the total pressure force on the window. The pressure 
exerted on the inside surface of the window was assumed to be freestream 
atmospheric. The dimensions of each rectangular area were determined by 
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equaling dividing the vertical and horizontal distances between each of 
the pressure ports located in the optical window area. 



Test Case # 


Window Force 


Force/Q-bar 


Test Case # 


Window Force 


Force/Q-bar 


1 


-51.3 


-1.51 


10 


-52.2 


-1.54 


2 


-40.4 


-1.19 


11 


-51.1 


-1.51 


3 


-40.4 


-1.19 


12 


-50.5 


-1.49 


4 


-37 


-1.09 


13 


-53.2 


-1.57 


5 


-48.7 


-1.44 


14 


-45.2 


-1.33 


6 


-41.1 


-1.21 


15 


-43.8 


-1.29 


7 


-44.9 


-1.32 


16 


-40.3 


-1.19 


8 


-46.2 


-1.36 


17 


-43.2 


-1.27 


9 


-50.7 


-1.50 


18 


-46.9 


-1.38 



Note: test case #18 conducted with antenna installed 



Table 3. Summary of Window Force for the 18 Test Cases 

Table 3 depicts the total pressure force exerted on the window 
(measured in lbf) for each of the 18 test cases. The window force varies 
from 37 lbf suction for test case 4 and to about 53 lbf suction for test case 
13. These total force figures were determined by extrapolating to full- 
scale. Also provided in Table 3 are the ratios of window force to Q-bar 
for each test case, so estimates at other freestream dynamic pressures can 
easily be made. 

C. CENTERLINE LONGITUDINAL PRESSURE VARIATION 

Sandia requested that pressure measurements also be taken on the 
upper fuselage section of the wind tunnel model. This information will 
assist them in their effort to properly position an external NACA vent to 
provide cooling for the onboard instrumentation. To accomplishment this, 
25 pressure ports were installed along the centerline longitudinal axis of 
the upper fuselage section at two-inch intervals, as measured along the 
surface. One wind tunnel test run was conducted at zero a/zero P with AP 
= 20.5 psf. 
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Figure 33 illustrates how Cp varies along the upper fuselage of the 
full-scale airframe. The shape of this curve is similar to that of a typical 
airfoil section. Stagnation Cp is observed near the leading edge and peak 
suction occurs approximately ten inches measured longitudinally aft of 
the leading edge. This information can assist Sandia in optimizing the 
placement of an external cooling vent. Ideally, this air exhaust vent 
should be placed in the forward fuselage surface section where suction is 
maximized. This results in maximum cooling airflow. However, even aft 
of the suction peak, the pressure remains below its freestream value, 
providing sources for external venting. 

Full Scale Upper Fuselage Cp vs Longitudinal Position 




Figure 33. Upper Fuselage Cp vs Longitudinal Position 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The surface pressures on the proposed aerodynamic fairing for the 
Altus II UAV have been successfully characterized. Using a half-scale 
model, a wind tunnel experiment conducted at one-half the Reynolds 
number experienced in full-scale flight has estimated the surface pressure 
coefficients. Two wind tunnel test runs were conducted initially to 
determine that this Reynolds number difference would be expected to have 
a minimal impact on data accuracy. 

Seventeen wind tunnel test runs, conducted at various a and (3 
combinations, illustrated the effect of small attitude changes on the 
fairing’s surface pressures. Of particular interest was the fact that the 
pressure contours over the fairing’s optical window remained relatively 
constant for these test runs. After measuring these window pressures, 
integration was performed to determine the estimated total force on the 
window area. Loads extrapolated to full scale varied from 37 to 53 lbf 
suction for a reference flight condition of 100 KTAS at sea level. An 
eighteenth test run determined that an antenna that was not available for 
the initial runs did not significantly affect the measured data. 

Lastly, pressure ports were installed on the upper fuselage of the 
model to take pressure measurements along its centerline longitudinal 
axis. These measurements illustrated the position of peak suction and 
characterized the flow aft of this peak. This information may be used to 
optimize the installation of external exhaust vents used for payload 
cooling. 

B. RECOMMENDATIONS 

This study provides baseline information for Sandia’s installation of 
the full-scale version of the proposed aerodynamic fairing onto the Altus 



39 



II UAV. The window forces that were calculated may be used to 
determine the required glass thickness needed to minimize weight, but 
also to minimize distortion. The pressure measurements on the upper 
fuselage may prove useful in optimizing the placement of external cooling 
vents. 

For the Naval Postgraduate School, a more permanent test 
equipment set-up should be acquired for future AWT testing. Virtually all 
of the equipment used for this thesis research was borrowed from other 
research areas. The AWT is an excellent research tool for low speed 
testing, but it needs dedicated test equipment to realize its full potential. 
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APPENDIX A. EXPERIMENTAL PROCEDURE CHECKLIST 



This section contains the checklist used to operate the wind tunnel 

and run each of the 18 test cases. A single test case took approximately 

40 minutes to conduct. The following checklist to conduct the wind tunnel 

experiment is presented below: 

1. Ensure the Solenoid Controller, SDIU, Signal Conditioner, PC 486 
and AWT Control Panel are turned ON. 

2. Set model AOA and Sideslip. 

3. Ensure the tunnel and test-section are clear of foreign objects and 
all access doors properly closed. 

4. Calibrate the SDIU/Scanivalve with the stand-alone H 2 O 
manometer. Using the Signal Conditioner, set the zero and span 
voltages which correspond to zero and 30 cm H 2 0 pressure. Zero 
pressure should read zero volts and 30 cm H20 pressure should read 
about negative 1.8 volts. 

5. Select port #45 on the SDIU to ensure sufficient overlap of 
collected data. Select remote operation mode on the SDIU. 

6. On the PC 486, NI Labview program starts automatically after 
bypassing the User Login. Select FILE - OPEN - “Vi. lib” - “SDIU 
4-20-99” - “SVWorking”. [Ref. 5] 

7. Once the “SVWorking” VI is activated, enter the file path for the 
3.5-inch floppy disk output data file. Enter “10” in the time delay 
field. Select ON for “Append to File”. 

8. At the AWT Control Panel, start the wind tunnel motor. 

9. After the AWT motor has stabilized, slowly advance the fan pitch 
angle and set 10 cm H 2 0 indicated dynamic pressure on the AWT 
H 2 0 manometer. 

10. Activate the Iterative run option on the “SVWorking” VI. 
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11. After “SVWorking” has measured and recorded the data from all 48 
ports, stop the VI. Change the pressure manifold on the Scanivalve 
to the next one. Select local mode on the SDIU, reselect port #46 
and return the SDIU to remote mode. Restart the VI to begin 
measuring and recording pressure data again. Repeat this procedure 
for each of the three pressure manifolds, for each test case. 

12. Use Notepad or another similar ASCII text editor to verify the data 
recording on the 3.5-inch floppy. Remove overlapping data 
elements and transfer the edited data into a Microsoft Excel 
spreadsheet program for further data reduction. 
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APPENDIX B. MICROSOFT EXCEL SPREADSHEETS 



This section includes the Microsoft Excel spreadsheets used for 
each of the 18 test cases. Also included is an expanded discussion of how 
Scanivalve/SDIU voltage measurements are converted to pressure 
measurements. The last two spreadsheets illustrate how the pressure 
forces on the optical window are integrated across its area to determine 
the total exerted force. 

For each test case, SDIU voltage readings are imported into the 
left-most columns of each Excel worksheet. Next, the voltages were 
converted to pressure measured in pounds/foot 2 using the following 
relations: 



Pressure (psf) = 



30 cmH 2 0 2116.2/75/ 



1033.68 cm// 2 0 



-1.8 Volts 

The calibration voltage of -1.8 volts corresponding to 30 cm H 2 0 on the 
stand-alone manometer varied slightly for each test run, but the basic 
conversion equation remained the same for all the runs. Recall that the 
pressure measured in Scanivalve port #1 was subtracted from all of the 
other 47 ports for all 18 test cases. This was the zero reference correction 
discussed in Chapter III Section B. 

In the last column set, the Cp is determined with the tunnel 
calibration and blockage correction applied according to the following 
equation: 



c,= 



P-Pj, 



[.9883 1 + .00 1 02 x (qbar )] x qbar 



V, 



, where Qbar = p^-p^ and 



V, — Vj 1 Model Frontal Area 

£-— = — X- 



So, therefore — = e + l. 

v 1 



Vj 4 Test Section Area 

The spreadsheets are presented in the following pages. 
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1.34E-01 2.38E-01 2.32E-01 
1.65E-01 2.38E-01 2 26E-01 
1.10E-01 2.35E-01 2 35E-01 
1.56E-01 2.53E-01 2.44E-01 
2.08E-01 2.41E-01 2.38E-01 
2.26E-01 223E-01 2.35E-01 
2.44E-01 2.17E-01 2.32E-01 
3.24E-01 2.14E-01 2.32E-01 
3.24E-01 2.04E-01 2.29E-01 
3.51E-01 2.1 IE-01 2.41E-01 
3.63E-01 2. 11 E-01 2.35E-01 
4 24E-01 2.17E-01 2.44E-01 
4.09E-01 2.29E-01 2.41E-01 
4.15E-01 2.23E-01 2.41E-01 
3.42E-01 2. 11 E-01 2.44E-01 
3.42E-01 2.08E-01 2.50E-01 
3.48E-01 2.1 IE-01 250E-01 
3.48E-01 2. 11 E-01 2.53E-01 

3.60E-01 2.08E-01 241E-01 
2.90E-01 2.08E-01 266E-01 

2 84E-01 2.08E-01 2.59E-01 
2.81E-01 2.14E-01 2.41E-01 
2.8 IE-01 2.17E-01 2.32E-01 
2.78E-01 2 1 IE-01 2 29E-01 
2.87E-01 2.14E-01 2.20E-01 
3.50E-01 2. 11 E-01 2.17E-01 
2.78E-01 2.14E-01 220E-01 
2.41E-01 2.17E-01 2.17E-01 
2.44E-01 2. 11 E-01 2 20E-01 
2.47E-01 2.14E-01 2.17E-01 

2.44E-01 2.20E-01 2 29E-01 
2.41 E-01 2.26E-01 2.14E-01 
2.44E-01 2.17E-01 2.14E-01 
2.53E-01 2.17E-01 2.47E-01 
2.81 E-01 2.23E-01 2.41E-01 

223E-01 

2.17E-01 

2.23E-01 

2 35E-01 


O.OOE+OO O.OOE+OO O.OOE+OO 

-5.74E-01 -5.62E-01 -5.83E-01 
-9.46E-02 1.59E-01 1.80E-01 

-1 46E-01 1 56E-01 1.68E-01 
-1.40E-01 1.37E-01 1.83E-01 

-6 41E-02 1.43E-01 165E-01 
-3.66E-02 1.59E-01 1.98E-01 

-2.14E-02 1.50E-01 1.77E-01 
O.OOE+OO 1.43E-01 1.68E-01 

1.83E-02 1.43E-01 1.68E-01 

5.19E-02 1.50E-01 1.83E-01 

7.63E-02 1.43E-01 2.04E-01 
1.25E-01 1.31 E-01 1.98E-01 
1.77E-01 1.37E-01 1.98E-01 
1.98E-01 1.34E-01 1.92E-01 

2.08E-01 1.37E-01 1 83E-01 
2.26E 01 1.37E-01 1.95E-01 

2.96E-01 1.34E-01 2.04E-01 

2.99E-01 1.25E-01 2.01E-01 

3. 11 E-01 1.34E-01 2.04E-01 

2.44E-01 1 40E-01 2.08E-01 

2 56E-01 1.28E-01 2.23E-01 

2.72E-01 1.40E-01 2.17E-01 

2.75E-01 137E-01 2.23E-01 

2.69E-01 1.25E-01 2.04E-01 

217E-01 1.28E-01 2.66E-01 

2.14E-01 1.13E-01 2.41 E-01 

2. 11 E-01 1.19E-01 2.20E-01 

2.17E-01 1.40E-01 2.17E-01 

2.08E-01 1.40E-01 2.23E-01 

2.14E-01 131 E-01 2.01 E-01 
2.17E-01 1.43E-01 1.92E-01 

1.92E-01 1.46E-01 2 20E-01 

1.74E-01 1.40E-01 0.2690 

1.80E-01 1.31 E-01 0.2297 

1.83E-01 1.53E-01 0.2493 

1.83E-01 1 43E-01 0.2802 

1.80E-01 1.56E-01 0.1933 

1.74E-01 1.56E-01 1.71E-01 
1.77E-01 1.62E-01 2.29E-01 

1.86E-01 1.62E-01 1.79E-01 

1.68E-01 
1.62E-01 
1 59E-01 
1.56E-01 
168E-01 
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6.71E412 8.85E-02 7.63E<>2 

•531 E-01 -5.31 e-01 -5 28E-01 
6.71E-02 217E-01 2 44E-01 
1.63E-02 2.41 E-01 2.53E-01 
O.OOE+OO 266E-01 2.66E-01 
7.63E-02 2.35E-01 2.26E-01 

1.43E-01 2.78E-01 2.99E-01 
1.59E-01 2.72E-01 2 93E-01 
2.04E-01 284E-01 2.93E-01 
1.86E-01 3.14E-01 3.11E-01 

2 20E-01 2.87E-01 3.14E-01 
2.29E-01 2.69E-01 3.20E-01 

3 11E-01 2.69E-01 3.1 IE-01 

3 33E-01 2.50E-01 3.05E-01 
3.72E-01 2.S6E-01 3.08E-01 
3.75E-01 2.47E-01 3.08E-01 

4 46E-01 2.59E-01 3.27E-01 
4.58E-01 2.93E-01 3.17E-01 
4.5SE-01 2.78E-01 3.24E-01 
4.15E-01 2.62E-01 3.24E-01 
4.03E-01 2.62E-01 3.27E-01 
4.06E-01 2.56E-01 3.36E-01 

4 21 E-01 2 44E-01 3.30E-01 
3.69E-01 2 47E-01 3.75E-01 

3.54E-01 2.50E-01 3.60E-01 
3.45E-01 2.62E-01 3.39E431 
3.42E-01 2.75E-01 3.39E-01 
3.48E-01 2 66E-01 3.42E-01 
3.51 E-01 2.62E-01 3.24E-01 
3.69E-01 2.89E-01 2.93E4J1 
3 27E-01 2.72E-01 3 24E-01 
3.1 IE-01 2.59E-01 3.48E-01 
3.14E-01 2.59E-01 2.72E-01 
3.1 IE-01 2.59E-01 2.81E-01 
3.14E-01 2.84E-01 3.39E-01 
3.05E-01 2 81 E-01 '2 26E-01 
3.11E-01 2.87E-01 2 26E-01 
3.14E-01 2.72E-01 2 87E-01 
3.4SE-01 2.84E-01 2.17E-01 
9.16E-02 2.78E-01 9.16E-02 
7.93E-02 2.78E-01 8 24 E-02 
7.32E-02 2.75E-01 8 24 E-02 
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5.19E-02 4.88E-02 5.19E-02 
-&WE-01 -5.74 E-01 -5.28 E-01 
3 36E-02 2 29E-01 2 53E-01 
0 OOE+OO 2 35E-01 2.78E-01 
-2 44E-02 2 56E-01 2 78 E-01 
2 44E-02 2 29E-01 2 62 E-01 

1.13E-01 2.56E-01 2.75E-01 
1.46E-01 2 56E-01 2 78E-01 
1.65E-01 2.56E-01 2 78E-01 
1.65E-01 2.89E-01 2.93E-01 
2 04E-01 2 62E-01 2.96E-01 
2 29E-0 1 2 44E-01 2 90 E-01 

2 75E-01 2 41E-0 1 2 96 E-01 

3 39E-01 2.44E-01 2 87E-01 

3.66E-01 2.41E-01 2.96E-01 

4 OOE-OI 2 35E-01 2 84 E-01 

4.73E-01 2 50E-01 3.08E-01 
4 76E-01 2.38E-01 3.08E-01 

4.18E-01 2 32E-01 3.05E-01 

4.12E-01 2.41E-01 3.30E-01 
4.09E-01 2.23E-01 3.14E-01 
4 21E-01 2.26E-01 3.1 IE-01 
3.72E-01 2 2OE-01 3.S7E-01 
3.51 E-01 2.29E-01 3 36E-01 

3 54E-01 2 41E-01 3.02E-01 
3 54E-01 2 38E-01 3.08E-01 
3.54E-01 2 38E-01 3.02E-O1 
3 60E-01 2 47E-01 2.72E-01 
3.27E-01 2SOE-01 3.0SE-01 

2 99E-01 2 41E-01 3.17E-01 
3.02E-01 2 4 IE-01 2.75E-01 
3.1 IE-01 2 41E-01 2.93E-01 
3.14E-01 2.44E-01 3.27E-01 
3.05E-01 2.50E-01 2 29E-01 
3.14E-01 2.53E-01 2.50E-01 
2.99E-01 2S3E-01 3 20E-01 

3 24E-01 2 56E-01 2.38E-01 

7.63E-02 2.56E-01 8.71 E-02 
7.02E-02 2 56E-01 6.71E-02 

8.71 E-02 2 66 E-01 8 7,E*2 
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•6.26E-01 -6,07E-01 .8.136-01 
O.OOE+OO 2.59E-01 1.98E-01 

2.756-02 2.53E-01 2.786-01 

4.58E-02 2.50E-01 2 32E-01 
8.246-02 2.6264)1 2596-01 
1.71 E-01 2.86E-01 262E-01 
2-17E-01 2.59E-01 2.386-01 
1.4864)1 2.53E-01 2.50E-01 
2.016-01 2.66E-01 2.72E-01 
2.53E-0 1 2.99E-01 2.7564)1 
2.75E-01 2 56E-01 2.81 E4>1 
2.78E-01 2.47E-01 2.75E-01 
4.08E4J1 2.47E4J1 2.5064)1 
4.18E-01 2.44E4)1 2.5364)1 
4.06E4)1 2.41E-01 2.62E4)1 
4.55E-01 2.29E4)1 2.59E-01 
5.0464)1 2.35E4M 2.7864)1 
4.91E4>1 2.50E-01 2.6964)1 
4.76E4I1 2.59E-01 2.69E-01 
4.276-01 2.44E4)1 2.6964)1 
4.1864)1 2.4164)1 2.75E-01 
4.1264)1 2.38E-01 3.05E-01 
4.1564)1 2.20E-O1 2.846-01 
4.0964)1 2.28E4M 2626-01 
3.78E-01 2.14E-01 3.1764)1 
3.426-01 2.17E-01 3.1764)1 
3.51 E-01 2.23E-01 2.62E-01 
3.33E-01 2.58E-01 2.9364)1 
3.57E-01 2.50E-01 2.5664)1 
3.48E4)1 2.356-01 2.4764)1 

4.06E-01 *3564)1 27864)1 
3.336-01 2.38E-01 2.596-01 

3.3064)1 2.3564)1 3.1164)1 
3.11E4)1 2.35E-01 2.7564)1 
2.96E-01 2.59E4)1 2.S9E-01 
3.1 1 E-01 2.59E-01 2.8764)1 
2.996-01 2.50E-01 2.906-01 
3.14E-01 2.446-01 262E-01 
8.85E-02 2.476-01 7.93E-02 
7.32E-02 2.41E-01 8.41E4I2 
7.02E-02 2.41E-01 8.10E412 
7.02E-02 2.38E-01 S.80E-02 


0.00fe*O0 4.^ 

-2,146-02 1.536-02 2.7564)2 
•6.0464)1 -8.82E-01 -5.49E-01 
-4 276-02 2.58E-01 2.38E-01 

3.97E-02 2.386-01 2.44E-01 

1. 166-01 2 47E-01 2 44E-01 
1.65E-01 2.476-01 2.4764)1 
1.196-01 2.386-01 2.47E4>1 
1.71 E-01 2.S0E-01 2.44E-01 

2 23E-01 2.596-01 2.62E-01 

2.506-01 2.356-01 2.69E-01 
2.566-01 2.23E-01 2.59E-01 
3.36E-01 2.236-01 2 44 E4)1 

3.576-01 2.17E4I1 2 53E-01 
3.7264)1 2.17E-01 2.47E4)1 
4.0064)1 2.14E-01 2.5664)1 
4.55E4)1 2 14E-01 2.5664)1 

4.49E-01 2.38E-01 2.53E-01 

3.78E-01 2.14E-01 2.686-01 

3.7564) 1 2.17E-01 2.7564)1 
3.88E-01 2.1 IE-01 2.78E-01 
3.91E-01 2.04E-01 2 476-01 
3.5764)1 2.1 IE-01 2.9964)1 
3.45E-01 2.04E-01 2.93E-01 
3.146-01 2.1 IE-01 2.53E-01 
3.24E-01 2.356-01 2.66E4)1 
3.30E-01 2.17E-01 2.50E4>1 
3.33E-01 2.146-01 2.26E-01 

3.396-01 2.236-01 2.47E-01 
293E-01 2.206-01 2.32E-01 
2.8464)1 2.14E-01 2.5664)1 
2.906-01 2.17E-01 2.S3E-01 

2.7564) 1 2.17E-01 2.4764)1 
2.84E-01 2.446-01 2.296-01 
2.816-01 2.29E-01 2.476-01 

2.906-01 2.26E-01 2.53E-01 
6.10E-02 2.29E-01 7.32E-02 
4.276-02 2.29E-01 5.80E-02 

3 66E-02 2.23E-01 5.80E-02 

3.666-02 2.23E-01 5.80E-02 
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-5.22E-01 -5.25E-0I >5.31 £-01 
-1.68E-0t 1.50E-01 1.59E-01 

-1.31E-01 1.53E-01 1.77E-01 

-1.07E-01 1.62E-01 1.63E-01 

-9.77E-02 1.31E-01 1.56E-01 

-641E-02 1.65E-01 1.95E-01 

1.53E-02 1.86E-01 2.04E-01 
4S8E-02 1.86E-01 2 23E-01 
8.55E-02 1.77E-01 2.1 IE-01 
1.07E-01 1.59E-01 2.14E-01 

1.07E-01 1.56E-01 2 08E-01 

1.86E-01 1.59E-01 2.1 IE-01 

2.04E-01 1.56E-01 2 04E-01 

2.20E-01 1.59E-01 2.14E-01 

2.32E-01 1.50E-01 2.11E-01 

2 90E-01 1.62E-01 2.14E-01 

3 02E-01 1.74E-01 2.1 IE-01 

3.1 IE-01 1.56E-01 2.20E-01 

2 87E-01 1.40E-01 2 20E-01 

2B1E-01 146E-01 2.32E-01 

2.84E-01 1.53E-01 2 38E-01 

2.84 E-01 1.53E-01 2.29E-01 

2 87E-01 1 53E-01 2.17E-01 

2 86E-01 1.48E-01 2.56E-01 

2 38E-01 156E-01 2.47E-01 

2.38E-01 1.65E-01 2.35E-01 

2 38E-01 1.53E-01 2.26E-01 

2.38E-01 1.50E-01 2.20E-01 

244E-01 1 S3E-01 2.14E-01 

2.44 E-01 1.50E-01 2.1 IE-01 

2 38 E-01 1.S9E-01 2.29E-01 

2.11E-01 1.62E-01 2.29E-01 

1.98E-01 1.62E-01 2 23 E-01 

2.08 E-01 1.71 E-01 2.20E-01 

2.11E-01 1 86E-01 2.3SE-01 

2.11E-01 1.74E-01 1.98E-01 

2.17E-01 1 80E-01 1.92E-01 

2.32E-01 1 83E-01 1.71 E-01 

1.89E-01 
1.95E-01 
1.98E-01 
2.04E-01 
2 08E-01 


-622E-01 -6.40E-01 -6.49E-01 
-1.16E-01 1.31E-01 1.59E-01 

-1 40E-01 1 S3E-01 1 83E-01 

-1.22E-01 1.56E-01 1.86E-01 

-9.77E-02 1.13E-01 1.43E-01 

-5.19E-02 1.59E-01 1 95E-01 

0 OOE+OO 1.77E-01 2.01E-01 
2 44E-02 1.77E-01 2.04E-01 
2.14E-02 1 80E-01 2 01E-01 

9.16E-02 1 89E-01 2.11E-01 
1.19E-01 183E-01 2.11E-01 

1.34E-01 183E-01 2.08E-01 

2.08E-01 1.92E-01 2.1 IE-01 

2.32E-01 1.92E-01 2. 11 E-01 

2 41 E-01 1.89E-01 2.14E-01 

2.50E-01 180E-01 2.1 IE-01 

3.1 IE-01 1.9SE-01 2.14E-01 

3 24E-01 1.89E-01 2.17E-01 

3 02E-01 1.77E-01 2.17E-01 

2 96E-01 1.83E-01 229E-01 

2.62E-01 1.86E-01 2.44E-01 

2.B2E-01 1.89E-01 2.3BE-01 

2.S9E-01 1.92E-01 2.38E-01 

2 66E-01 1.89E-01 2 29E-01 

2.59E-01 1.92E-01 2.17E-01 

2.69E-01 1.92E-01 2.17E-01 

2 47E-01 1.95E-01 2.3SE-01 

2.29E-01 1.95E-01 2.23E-01 

2 44E-01 2.14E-01 2.66E-01 

235E-01 2.04E-01 1.92E-01 

2 32E-01 2.14E-01 1.80E-01 

2.32E-01 2.08E-01 2.41E-01 
244E-01 2 20E-01 189E-01 

2.20E-01 

2.23E-01 
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1 6.10E-03 3.36E-02 1.83E-02 

2 O.OOE+OO O.OOE+OO O.OOE+OO 

3 -5.58E-01 -5.62E-01 -S.58E-01 
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APPENDIX C. MATLAB GRAPHICS PROGRAMS 



This section contains the MATLAB programs that transformed the 
Microsoft Excel spreadsheet data into the color pressure contour charts. 
These programs were written by Professor Howard to maximize the 
visualization of the collected data. A summary of these MATLAB 
programs follows: 



1. ‘conplot’ plots fairing, window & color pressure contours 

2. ‘dots’ plots fairing, window & port positions 

3. ‘fairing_plot’ data file for fairing plot 

4. ‘window’ data file for window plot 

5. ‘data824’ data file of Cp’s for ‘conplot’ 



The data file (e.g. ‘data824’) changed periodically as data reduction 
matured. The ‘824’ portion of the filename referred to the date 24 August 
1999, the date of the latest reduced data. The MATLAB codes for the 
‘conplot’ and ‘dots’ programs are presented in the following two pages. 
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MATLAB PROGRAM ‘CONPLOT 



load data824.m -ASCII 
a = data824; 
x = a(:, 1 ); 
y = a(:,2); 
for i = 1:18 
z = a(:,i+2); 

xlin = linspace(min(x),max(x), 100); 
ylin = linspace(min(y),max(y),100); 
[X,Y] = meshgrid(xlin,ylin); 

Z = griddata(x,y,z,X,Y, ’cubic’); 

contour(X,Y,Z,15) 

caxis([-0.6 0.1]) 

colorbar 

%axis off 

% 

hold on 

load fairing_plot.m -ASCII 
b = fairing_plot; 
x2 = b(:,l)/2; 
y2 = b(:,2)/2; 
plot(x2,y2,’k’) 

xlabeI(’Model Dimension, inches’) 

ylabel(’Model Dimension, inches’) 

title('Pressure Coefficients’) 

load window. m -ASCII 

a3 = window; 

x3 = a3(:,l ); 

y3 = a3(:,2); 

x3 = x3/2; 

y3 = y 3/2; 

plot(x3,y3,’k’) 

%axis U-15 15 -30 25J) 

%axis equal 
axis ([-30 30 -30 30]) 
hold off 
% 

%axis off 

figure 

end 



62 



MATLAB PROGRAM <DOTS 



load data824.m -ASCII 
a = data824; 
x = a(:,l); 
y = a(:,2); 
plot(x,y,V) 
axis ([-30 30 -30 30]) 
axis equal 
% 

hold on 

load fairing_plot.m -ASCII 
b = fairing_plot; 
x2 = b(:,l)/2; 
y2 = b(:,2)/2; 
plot(x2,y2,’k') 

xlabel(’Model Dimension, inches’) 

ylabel(’Model Dimension, inches’) 

title(’Pressure Port Locations') 

load window. m -ASCII 

a3 = window; 

x3 = a3(:,l); 

y3 = a3(:,2); 

x3 = x3/2; 

y3 = y 3/2 ; 

pIot(x3,y3,’k’) 

%%a.xis equal 
hold off 
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